Thin Shell Dynamics and Equations of State 
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A relation betweeen stress-energy and motion is derived for accelerated Israel 
layers. The relation, for layers between two Schwarzschild manifolds, generalizes the 
equation of state for geodesic collapse. A set of linked layers is discussed. 
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I. INTRODUCTION 

Equations of state and boundary matching are two important tools for developing exact 
solutions of Einstein's field equations. As models have become more physical, the require- 
ments of strict boundary matching have been relaxed by joining two exact solutions across 
a boundary layer, matching their metrics on the layer but allowing jumps in the derivative 
structure. The Israel junction conditions [l] are often used to make the broader matching 
and are widely applied because they provide a simple dynamic boundary description for a 
variety of scenarios ranging from thin shell descriptions , to shell applications like 
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ce dilatons 
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141 ]. gravastars 



and extensions of 
20j or Gauss-Bonnet gravity |2ll. l22j. Since their introduc- 



have played an increasingly important role in gravitational physics. 
Barrabes and Israel [23j] began their paper with a description of the Israel layer as a ther- 
modynamic phase boundary, but the initial applications of Israel layers considered metric 
matchin g in dynamic collapse processes involving dust shells, null shells, and cosmic string 
loops [l|, 



Poisson 



24J has summarized some of the early seminal work 



1. 



23j, |25|, [26]. 



As our knowledge of the variety of astrophysical objects and their dynamic processes has 
expanded, Israel layers have become physically interesting in their own right and the ques- 
tions that have been investigated for large scale three-dimensional mass distributions, are 



now being asked about Israel layers 27J, 



28 



29 



30. 
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In this paper we investigate the relation between a layer's equation of state and its motion. 
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An equation of state for a layer dropping from rest at infinity on an exterior Schwarzschild 



geodesic has been formulated |32| , and can be generalized to include accelerated motions. We 
develop the extension and apply it to static and dynamic examples. The motion input to the 
geodesic extension is a single function. In the next section we briefly review the thin shell 
formalism used in the rest of the paper and give the geodesic extension. The applications of 
the extension are in Section III. Many of the applications use a layer with a linear equation 
of state, P = ao (L-layer). The description of an L-layer has been included in Appendix 
B. The accelerated extension singles out several special 'a! values, a = (—1/2, —1/4, 0). We 
show that the a = —1/4 L-layer is simply related to a geodesic layer. A set of linked L-layers 
is described and the other two special 'a' values are shown to be boundaries for the linked 
layers. 

II. THIN SHELL DESCRIPTION 

We consider an Israel layer as a thin shell S between two Schwarzschild manifolds with 
exterior mass parameter m and interior mass parameter M. The spacetime consists of the 
two manifolds which join across surface E (+/- denotes exterior /interior). 

g±dx a dx b = -{f±)dt 2 + (l/f±)dr 2 + r 2 dtt 2 , (la) 
/+ = l-2m /R, f. = l-2M/R, (lb) 
dsl = -dr 2 + R 2 (r)dn 2 . (lc) 

The layer is tracked by two observers comoving with the layer. The observers use r = 
R(t) and T± = T±(t) to describe the layers. With this parametrization and with overdots 
denoting d/dr, the observers' velocities and associated normal vectors are 

Ui={T±,R,Q,0), i = t,r,#,<p (2) 
n i± = (-R,f ±,0,0). (3) 

Velocity normalization g^U±U± = —1 implies (f±) 2 T 2 = R 2 + f±. The jump in f±, T±, is 
an important input to layer motion, and we define 

A = A+ - A_ (4) 



A± = /±T ± = V^ 2 + /±- (5) 
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Matter Content 



The layer has stress-energy content 

S\ := aU% + Ptfj + 17%) 

with density a and stress P ( — P is tension). The stress-energy content o 
on the bounding metrics and the layer motion. From the Israel conditions 

4trj = -A/R, 
8nP = A/R + A/R, 



(6) 

? layers depends 
24J we have 



with the layer mass defined by 



mi = 4ttR 2 it. 



(7) 
(8) 

(9) 



Layer Motion 

Motion parameters R 2 and R are both important inputs to the stress-energy structure of 
the layer. 



' mo-M V | mp+M 



R 



A 

7 



i?A J R 
(m -M) 2 ] m + M (m — M) 2 



A 3 R 2 



2R 2 



A 2 R 3 



(10) 
(11) 



The radial components of the 4-accelerations of observers comoving with the layer (Appendix 
A) are 

ir ± = R+(l-f±)/(2R). (12) 
The radial accelerations are related to a stress-energy sum 



4vr(a + 2P) 



A 
R 



u r 



(13) 



T + R 2 -2m /R + R 2 -2M/R 

These relations provide insight about the role of jumps in the observer accelerations and 
velocities in defining the stress-energy. It is clear that the stress-energy structure of the 
layer and its motion are related so that different assumptions about the motion will produce 
different equations of state or inversely, that an imposed equation of state will determine 
the layer motions. 
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The Geodesic Generalization 

A layer starting from rest at infinity and collapsing along a geodesic in the exterior 
spacetime has motion parameters 



R g = -^/2m /R, (14a) 
R g = -mo/R 2 , (14b) 
i? 9 (r) = (2m ) 1 / 3 [(3/2)(c 1 -r)] 2 / 3 . (14c) 



The equation of state for the geodesic layer [32J is 



I'^W'-W'S' (15) 

This geodesic equation of state is the limit (7 — > 0, A — > 0) of a layer with exterior motion 
parameters 

where A is the radial 4-acceleration of the exterior co moving observer and 1 + 7 = A+. A 
and 7 are related by A = j/(2R). With these parameters, Eq. (fT3]) can be written as 

Ma + 2P) = * - A±B^m . (18) 
VTT^ y/l+j + 2(m -M)/R 

In this equation, the motion parameters, R, A, and 7 require several separate choices to set 
integration constants. The stress-energy-motion relation (Appendix C ) involving a single 
motion function F := ARj (1 + 7) is 

7r(x(m - M)(l + AP/af = (1 + 7 ) 3/2 (F - 2P/a) 2 (1 + F + 2P/a). (19) 

This is the generalization of the geodesic equation of state. For m = M there is no layer 
and the motion is described by RA = 0, 7 = const. In developing applications of this 
equation, one notes that two key inputs are either the ratio Pj a or the value of F. Choosing 
a static layer (R = 0, R = 0) sets the value of F and a general equation of state results. 
Choosing P/a selects the motion; for example, the equation of state Pjo = (—1/4, —1/2, 0) 
are all simplifying special values for the relation. In the next section we will examine all 
three values. We begin with dynamic layers and consider static layers as a second example. 
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III. APPLICATION TO NON-GEODESIC LAYER MOTION 

A. Dynamic Layers: P=-(l/4)<7 

P = ao = — (l/4)cr is singled out by the generalized stress-energy-motion relation. From 
Eq. ffT9l . for this equation of state, one has 

AR 



1 + 7 

This can be integrated for 1 + 7 giving 



-1/2. (20) 



l + 7 = f (21) 
with Cq an integration constant. Using the motion function, A is 

A = = R g . (22) 

2R 2 2m 9 V ' 

Layers that have a radial 4-acceleration linearly related to the geodesic R g value will have 
the P = —a /A equation of state. They are L- layers with tension. In the section on static 
layers, we will see that a = —1/4 is an excluded value for static L-layers. We will also 
see that a = —1/4 is an important boundary point for some of them. Note that there are 
moving L-layers with tension. 

B. Linked dynamic layers 



Because of its simplicity, the linear equation of state is often used in discussing thin shells 



33(. For L-layers (Appendix B) 

A = - Ca R- {1+2a) . 



For this A, the motion of the layer is described by Eq. fflO]) 

= clR-^a) ( mo - M f m + M _ 

4 c 2 a R- Aa R 1 ' 

There is an interesting symmetry in this equation. Over the positive 'a! range, a = a p , 
< a p < 1 Eq. (j2U|) describes the motion. For a = —a n , < a n < 1 the motion is described 
by 

= ^-2(1-200 (mq _ M) 2 mQ + M _ 

" ' 4 R 
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The substitutions 

a n = a p + 1/2, (25) 
c 2 n = 4(m -M) 2 /c 2 p , 

map these two motion equations into each other. The 'a' range over which a layer with 
tension is linked to a layer with pressure is 

l/2<a„<l (26) 
< a p < 1/2. 

For each negative a = —a n in this range, there is a layer with tension which has the same 
motion as a layer with pressure and positive a = a p . The linked layers have the same R and 
R. a p = and a n = 1/2 are the lowest 'a' values for linked pressure/tension shells and are 
two of the values which simplify Eq. (fl9l) . For these two values we have 

a = : 7r(m - M)a a=0 = (1 + 7 ) 3 / 2 F 2 (l + F) (27) 
a = -1/2: 7r(m - M)a a= _ 1/2 = -(1 + 7 ) 3 / 2 (F + 1) 2 (F) (28) 

The motion functions are 

r 2 

a = : F — 



-1/2: F 



R[cl - 2R(m - M)] 
2(m - M) 



R[2(m — M) — Rc 2 _ 1/2 ] ' 
and are identical under the linkage. The matter content of the two layers is different. From 
Eq.(|271) and Eq.flSg) we have 



O"a=0 



Va=-l/2 1 + F 

For m = M, we would expect no layer to exist and a = —1/2 describes Schwarzschild 
vacuum. 

A missing part of the 'a' range for linked pressure/tension shells is— 1/2 < a < 0. There 
are linked L-layers in this region but the linked layers both have tension. This linkage is 
centered around a = —1/4, one of the special values for the geodesic extension, describing 
layers whose radial 4- accelerations are linearly related to R g . Consider two negative 'a' values, 
a ni and a„ 2 in the ranges 

1/2 < a m < 1/4 (29) 
1/4 < a n2 < 0. 
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As might be expected from the pressure/tension linking relations, the relation for these 
layers are 

a ni = -a n2 + 1/2 (30) 

c 2 anl =4(m -M) 2 /cl 2 . (31) 

When a = —1/4, then a ni and a n2 coincide. 



C. Static layers 

Static layers with R = 0, R = have motion parameters 



7s = -2m /R s , (32) 
A s = m /R 2 s 



1 S ! 



and are, in some sense, the negatives of the geodesic layer with j s = —R 2 ,, A s = —R g . Using 
Eq.© and Eq.flTSJ) we have 

4P/a = -1 + (33) 

1 - 2M/R S - 4irR,<Ty/l - 2M/R S 

For M = this is the stress energy-radius relation for static layers given by Khourami and 
Mansouri |3j, with their stress and density related to P and a by P^m = 8nP, a^m — 8na 

P = ^ . (34) 

This relation is particularly interesting when compared with the classical van der Waals 
form for 3-dimensional fluids. 

nRT . n s 9 

with n the number of moles, volume V, gas constant R, and temperature T. a and b are 
constants. The denominator of the first term corrects for a minimum volume available to 
the fluid constituents. This can be attributed to a finite size of the particle constituents or 
to the existence of a repulsive core in constituent interactions. The second term accounts for 
an attractive long range attraction between constituents which reduces the stress. For low 
densities (n/V) the equation of state describes a perfect fluid, PV = nRT. 

The numerator of Eq. (l34l) could imply that, for irR s a << 1, an Israel layer is a first 
order polytrope. The current relativistic polytrope assumes a linear low density equation of 



8 

state [5]. The part of the classical van der Waals equation describing long range interaction 
is missing from Eq. fl34l) but the denominator suggests, as in the classical equation, there 
is either a minimum or a zero. The existence of a minimum value could be related to the 
existence of a repulsive core in the interaction potential between allowed layer constituents. 
Detailed models of this possibility will be discussed elsewhere. Equation ( 1341) also suggests 
that there are static shells with tension and pressure with the zero denominator related to 
the boundary between the two kinds of stress. At the zero value the Israel A can be found 
from the density 

AnR s a = 1, 
A = -1. 

R = is one of the static shell conditions. From the definition of A, Eq.(j4|), we have for an 
M = static shell 

A=Jl 2 '"" 



R 

We see that there are no static shells outside the horizon corresponding to A = —1. The 
boundary between the static linked layers is not static. 

The static shell boundaries can be explored using M = 0, L-shells as an example. For 
static L-shells the Israel radius value is (Appendix B) 

R -m {Aa+1)2 
Rs ~ m °4a(l + 2a) 

and the range for physical static radii is —1 < a < 1/2, < a < 1. There are linked static 
shells with tension and pressure, just as in the dynamic case. For each R s in the range 
< a < 1/2, there is an identical radius in the range — 1 < a < —1/2. For example, 
a = 1/2 and a = — 1 have the same static radius R s = |m , as suggested by the linking 
relation a n = a p + 1/2. The range for static shells excludes the points a = — 1/2 and a = 0. 
These points are the moving boundaries of the static layer region. The boundary between 
the pressure/tension linked M = static L-layers is a moving layer that is not the dust layer 
one might have expected, but the layer linked to dust. 

The M^O static relation is less easy to interpret in a van der Waals sense, since the size 
of the interior mass M becomes an important parameter. 

P ^ 2 | M ^/ 2 (35) 

y/l - 2M/R S - A-nR s a R s - 2M - AnR s a^/R s - 2M ' 
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For AkR s o « -i/l — 2M/R S the equation of state becomes approximately linear 

2M , . 

P « a (36) 

R S -2M K J 

and a mixed equation of state results as the effect of M increases. 

IV. CONCLUSION 

A general stress-energy-motion relation was derived for Israel layers between 
Schwarzschild manifolds. The equation generalizes the equation of state for layers drop- 
ping on exterior Schwarzschild geodesies. It was used to discuss the relation between motion 
and equation of state. The motion input is a single function of the exterior comoving ob- 
server acceleration and velocity. Using the relation, the motion of a layer with equation of 
state P = —a /A was shown to be linearly related to geodesic motion. 

A set of linked L-layers with a common motion was described over the parameter range 
— 1 < a < —1/2, < a < 1/2. In the linked range, each layer with pressure has a partner 
layer with tension. There is also a set of linked shells, both with tension in the range — 1/2 < 
a < 0. These layers coincide for P = — cr/4. There are layers in the range 1/2 < a < 1 but 
they are not linked to physical layers with tension. Positive 'a' values in this range are 
linked to negative values larger than 1. Because the motion of the linked shells is the same, 
their motion functions agree and the new relation can be used to compare densities and 
pressures. 

The stress-energy in the Israel formalism is described by two observers comoving with 



HQ , 



the existence of static layers with 



the layer. As pointed out by Ipser and Skivie, 
tension is related to the accelerations of the comoving observers. For a static layer, the 
two observers are hovering over the layer but must accelerate in order to remain static with 
respect to the layer. For L-layers, the radially projected 4-accelerations given in Eq. fflBl are 
related to the size of state parameter 'a'. 

n r+ U r + - n r Jj r _ = 4vra(l + 2a). (37) 

This can be used to interpret the relative sizes of U r needed for the two hovering observers. 
For the M = case, the interior observer is not accelerated at all and this equation describes 
whether the projected 4-acceleration needed by the exterior Schwarzschild observer points 
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inward or outward. For — 1/2 < a the observer needs to accelerate away from the layer, 
counteracting the gravitational attraction of the layer and for a < —1/2, the region where 
there are static shells with tension, the observer has to accelerate toward the layer. The 



repelling nature of layers with tension has also been discussed by Vilenkin [34j • 

The applications focused mainly on L-layers because the P/a structure of the geodesic 
extension makes this a simple example to develop with clarity. There are many interest- 
ing questions yet to be studied. Other equations of state easily could be investigated, for 
example, a dynamic first order polytrope with P = Ko 2 would have a density 

, (1+4 ^ ) 3. (1+7 p ( ^_ 2&) ^l+™\ 

Using this equation, the suggestion that low density static shells are first order polytropes 
could be explored for dynamic shells. The general stress-energy-motion equation has been 
developed for Schwarzschild but could be a useful tool in understanding layers bounding 
other metrics. The layer linkage discussed here depends on the radial structure of A in 
Schwarzschild L-layers. The idea of linkages for L-layers bounding other metrics in a variety 
of dimensions may have broad applications. 

APPENDIX A: OBSERVER ACCELERATION 

The layer is tracked by two observers comoving with the layer who agree on the layer 
metric. The radial 4-accelerations of these observers are computed from the 4-acceleration 

A* = U^aU 1 = U a d a U l + T^tTfA (Al) 

Using 

dR dRdt dRdr dR- OR- 

_|_ rp _|_ 

dr dt dr dr dr dt dr 



one finds 



A r = § + i^f](Ty-Q(l/f)}(Rf. (A2) 



The velocity normalization is 

ft 2 = R 2 /f + 1. 



For the exterior observer we have 



* r = R + l£- ( A3 ) 
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Thus, in general 

U r ± = R+(l-f±)/(2R). (A4) 

APPENDIX B: P = aa 

General Motions 

Because of its simplicity, the linear equation of state is often used with the field equations 
and it is frequenly used in discussing thin shells . For layers with P = aa, we have 

£ = -(l + 2.)§. (Bl) 

If 'a' is a constant we have 

A = -c a R- {1+2a) . (B2) 
AnP a = ac a R- 2{l+a \ 
4na a = c*R- 2 < 1+a l 
m La = c a R~ 2a . 

c a carries an 'a' index because, from a unit standpoint, it will have to vary with the value 
of 'a'. The motion of the layer is described by 



c 2 a + (m - M) 2 R Aa + m + M _ ^ (fi3) 



4jR 2(l+2a) c 2 R 

» c 2 a {l + 2a) 2a(m -M) 2 R^- 1 m + M 

R= ~ 4i?3 + 4a + -2 2iP~' (B4) 



'a 



The points R — provide an equation for c a . 



2Rl +4a [R -m -M± yJR* - 2^0(^0 + M) + 4m M] = < (B5) 

Not all values of 'a' will correspond to a static layer with both R — and i? = 0. Where 
there is a static layer, the constant, be evaluated in terms of the equation of state 

parameter 'a'. 
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Static layer 



The static points, R s , follow from R s = 0, R s = 0. This identifies Rq with R s . Using 
Eqs.( TB4|B5ll one finds 

(1 + 2a)4R 7 ^ - 8aim °- M » 2 + 2 ''"° + = p. (B6) 

The parameter values, a = 0, a = —1/2 have no solutions. The static radius in terms of c a 



is 



R i+4* = (1 + 2Q)C " . (B7) 

(rn + M)± y/(mo + M) 2 + 8a(l + 2a) (m - M) 2 

Using Eq.( lB5l) for c a , The static layer radius is 



(m + M)(l + 4a) 2 ± ^(toq + M) 2 (l + 4a) 4 - 32am M(l + 2a)(l + 4a) 2 

^ = 8a(l + 2a) " (B8) 

For M = the layer description is especially simple and we have 

(4a + 1) 2 , , 

R = m ±- J —. B9 

°4a(l + 2a) V ; 

2 _ ci+4a 4am o 



c: = i2 



l + 2a 



4- a ^ ±M!il±M. (BIO) 



4vrP a 



m (4a + l) 3 
1 16a 3 (l + 2a) 



m (4a + l) 3 
1 + 4a 



over the 'a' range 



-1 < a < -1/2 
< a < 1 

It is clear that there are static shells with tension as well as with pressure. 



APPENDIX C: DERIVATION OF GEODESIC EXTENSION 

The extension is developed in terms of A and 7, which are zero for exterior geodesic 
motion. 

A = R + m /R 2 . (Cl) 
7 = R 2 - 2m /R. (C2) 



The Israel formalism gives the stress and density in terms of a function 



8nP = A/R + A/R, 
Ana = -A/R, 



with 



A+ = Jl + R 2 - 2m /R = V 1 + 7- 



A_ = y/l + B?-2M/R = ^l + 1 + 2 m ° R M . 
A = A+ - A_ 

Calculating the derivatives with respect to r we have 

R(R + m /R 2 ) RA 



A, 



l + R 2 - 2m /R ^ 



A _ R{R + M/R 2 ) _ R(A + M/R 2 - m /R 2 



l + R 2 - 2M/R A - 
There are several useful relations. Using Eq. flC4j) and (IC5I) one finds 

Al - A 2 = 2- 



A 2 „M-m 



R 

and using this, the density of the layer is 



Ana 



(A + -A_)(A 2 + -A 2 



2(m - M) 

A pressure-density relation can also be found 



8ttPR + AiroR 
8nPR + AixoR 
8ttPR + AnoR 



RA , 
-AnaR + | 




RA(A_ - A+) 


M -mo 


A_A+ 


RA ' 


iL4(A_ - A+) 


A+ - A 2 


A_A+ 


2A_ ' 


RA{AnaR) (A. 


f + A_)AiraR 



A_ 



2A_ 



AR A + 



2P/a + 1/2 = — — + 



9 AR 

{APj a + 1)A_ = + A + . 



A_A + 2A_ 
IA. 

a" 
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A useful relation is 



P 

(4- + 1) 



A 
A 



+ 



2AR 
1+7 



+ 1 



(C9) 



with 





A_ 



2AR/(l + 7 ) + 2 + 4P/or 

(AP/a + 1) 
2AR/(l + 7 ) -4P/cr 
(AP/a + 1) ' 



(ClOb) 



(ClOa) 



Using these, substituting into Eq.( lC8j) . the general extension of the geodesic equation of 
state follows. 



For A = 7 = 0, this becomes the geodesic equation of state. 
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